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We investigate the ground state and the thermal entanglement in the two-qubit Ising model inter- 
acting with a site-dependent magnetic field. The degree of entanglement is measured by calculating 
the concurrence. For zero temperature and for certain direction of the applied magnetic field, the 
quantum phase transition observed under a uniform external magnetic field disappears once a very 
small non-uniformity is introduced. Furthermore, we have shown analytically and confirmed numer- 
ically that once the direction of one of the magnetic field is along the Ising axis then no entangled 
states can be produced, independently of the degree of non-uniformity of the magnetic fields on each 
site. 

PACS numbers: 03.67.-a,03.65.Ud,75.10. Jm 

I. INTRODUCTION 

The existence of entangled states of component systems and their unique properties have attracted a lot of atten- 
tion since the early days of quantum mechanics [1-5]. Entanglement has also recently been the subject of several 
investigations as it plays an important role in the topical areas of quantum computation and quantum information 
processing [6,7]. 

The Heisenberg magnetic spin system is one of the few physical systems that entanglement arise naturally. Nielsen 
[7] was the first to report results on entangled states utilized in a two-spin system. He calculated the entanglement of 
formation of a magnetic system described by an isotropic XXX Heisenberg Hamiltonian in an external magnetic field 
directed along the z-axis. Later, Arnesen et al. [8] systematically investigated the dependence of entanglement on 
temperature and on the applied external field in the same ID Heisenberg system. They showed that for ferromagnets 
the spins are always disentangled while entanglement is observed for antiferromagnets. In addition, they found that the 
degree of entanglement can be enhanced by increasing the external parameters (magnetic field and/or temperature). 
Then, Wang [9,10] studied the effect of anisotropy on thermal entanglement, working on the two-qubit quantum 
anisotropic Heisenberg XY model [9] and on the anisotropic Heisenberg XXZ model [10]. He also investigated the 
isotropic Heisenberg XX model with an external magnetic field applied along the z-axis [9]. Kamta and Starace [11] 
investigated the thermal entanglement of a two-qubit Heisenberg XY chain in the presence of an external magnetic 
field along the z-axis. They showed that by adjusting the magnetic field strength, entangled states are produced for 
any finite temperature. Sun et al. [12] extended later the work reported in Ref. [11] by introducing a non-uniform 
magnetic field. Comparing to the uniform field case, they showed that entanglement can be more effectively controlled 
via a non-uniform magnetic field. The full anisotropic XYZ Heisenberg spin two-qubit system in which a magnetic 
field is applied along the z-axis, was studied by Zhou et al. [13]. The enhancement of the entanglement for particular 
fixed magnetic field by increasing the z-component of the coupling coefficient between the neighboring spins, was 
their main finding. Finally, thermal entanglement in a two-qubit Ising model assuming an applied magnetic field in 
an arbitrary direction has been investigated by Gunlycke et al. [14]. 

Several, mainly numerical, works exist also on the study of pairwise thermal entanglement in the n-qubit Heisenberg 
spin chain. In short, the systems being studied are the XXZ three-qubit Heisenberg model with an applied magnetic 
field in the z— direction [15], the XYZ three-qubit Heisenberg model with an applied magnetic field along the z-axis 
[13], the XX four-qubit Heisenberg model with an applied magnetic field in the z-direction [16], the XX [17] and XXX 
and XXZ [18] four- and five-qubit Heisenberg model and the XX m-qubit (with m up to eleven) Heisenberg model 

It turns to be quite interesting to study the general case of different magnetic fields at each spin site. The control 
of the applied field at each spin separately is very useful in order to perform quantum computations [20]. Hence in 
the present theoretical analysis we investigate the ground state and the thermal entanglement in an Ising model with 
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a non-uniform and anisotropic external magnetic field, i.e. the magnitude and the direction of each magnetic field is 
different in each spin. 

This article is organized as follows: In the following section we present the details of the theoretical analysis based 
on the calculation of the concurrence [21,22] of the system. Then, in section III numerical results are presented 
and discussed for several cases of the parameters of the system (magnetic fields magnitudes and directions and 
temperature). In addition, a general analytical result is presented in the Appendix. There, we show that once one of 
the applied external magnetic fields points along the Ising direction, no matter what is the direction and magnitude 
of the other field the concurrence is always zero. This result is valid for zero and finite temperatures. Our results are 
summarized in section IV. 



II. THEORY 

The Hamiltonian studied in this work is given by 

H = J(o{o z 2 +<y^al) +B 1 cos6»ictJ + Bisin0iCT^ + B 2 cos 6» 2 <r| + B 2 sin 9 2 <7 2 > (!) 

where a a are the Pauli operators and J is the strength of the Ising interaction. Also, B\ and B 2 are the magnitudes of 
the external magnetic fields. We assume that each magnetic field has an arbitrary direction, defined by the angles Q\ 
and 9 2 between the field and the Ising direction. It is sufficient to consider that the magnetic field lies in a plane (xz) 
containing the Ising z-direction, because in three spatial dimensions, the Hamiltonian possesses rotational symmetry 
about the z-axis. 

A useful and convenient quantitative tool, which has been developed to study entanglement, is the entanglement 
of formation [4] . The entanglement of formation of a state of a composite system is proportional to the minimum 
number of Bell states, which must be shared between the components of the composite system. There is no general 
prescription for evaluation of the entanglement of formation in arbitrary systems. In this area Wootters [21] has 
described an explicit method for evaluating the entanglement of formation of an arbitrary state of a two-qubit system. 
His explicit formula is a generalization of the expression derived by him and Hill [22] for a special class of density 
matrices. 

For a pair of qubits the entanglement of formation E\ 2 is estimated from the expression E\ 2 = h (l + yjl — Cf 2 ) /2, 

where C\ 2 denotes the concurrence and h is the binary entropy function [21]. The concurrence is defined as C12 = 
max{Ai — A2 — A3 — A4, 0}, where the A's are the square roots of the eigenvalues in decreasing order of magnitude (i.e. 
Ai > A2 > A 3 > A4) of the spin-flipped density matrix operator R 12 = p{<J v <g> <J v )p * (<J y (g> <r v ). Here, p is the density 
matrix operator defined as p = ex-p(—H/ksT)/Z, where Z is the partition function (Z = tr{cxp(— H/ksT)}) and ks 
is the Boltzmann's factor. As concurrence is a monotonically increasing function of E\ 2 and both functions have values 
in the range to 1, we practically use C12 as a measure of the entanglement. Zero concurrence corresponds to an 
un-entangled pair of states and unity concurrence to a maximally entangled pair of states. This type of entanglement 
is usually called thermal entanglement as it is described by a temperature dependent density matrix operator. 

In order to proceed we need to find the eigenvalues (Ei) and eigenvectors of the Hamiltonian of the Ising 

system in the presence of a non-uniform external magnetic field, i.e. the Hamiltonian of Eq. (1). Once we have 
determined the eigenstates of the system, the density matrix operator can be written as 

p = Z- 1 Y,cM-E l /k B T)\V l )(* l \ . (2) 

i 

Then, the spin-flipped matrix operator is evaluated in a 4 x 4 matrix representation, in terms of the natural basis 
vectors {|00) , |01) , 1 10) , |11)}. In most cases, even if one obtains analytic expressions for the eigenvalues of the spin- 
flipped matrix operator, it is practically impossible to derive a simple analytic expression for the concurrence. The 
reason lies to the fact that the relative order of magnitude of the eigenvalues of R\ 2 , depends on the parameters 
involved. In general the concurrence can be evaluated numerically. For particular values of the parameters of the 
system, an analytic expression can be achieved (see in the Appendix). It is worth pointing out that the above 
analysis is greatly simplified for the zero temperature case, and then only the ground state is populated and hence 
p = l^o) (*o|; where is the index that denotes the ground state of the system. 

III. RESULTS AND DISCUSSION 

We begin our discussion with the results of the ground state of the system at zero absolute temperature (T = 0). 
In Fig. 1(a) we assume a uniform magnetic field (i.e. same magnitude and same direction of magnetic field in each 



spin) and present the results of the concurrence as a function of the strength of the applied fields. It is clear that the 
entanglement is highest for nearly vanishing magnetic fields and decreases with increasing the field's amplitude. In 
this figure we observe that for B\ = B2 = B — > + , the concurrence is unity, which shows the creation of maximally 
entangled states. Once we set a zero value for the magnetic field, the eigenstates are the same as those of the 
Ising model without magnetic field, i.e. the standard disentangled basis {|00) , 1 01) , 1 10) , |11)}, as the Hamiltonian 
is diagonal in this basis. Therefore, the concurrence obtains a zero value. Then, even for an infinitesimal increase 
of the magnetic field, the system goes from a non-entangled state to a maximally entangled state. This is a clear 
evidence of a quantum phase transition (QPT) [23]. In the case that the magnetic field is along the z-axis, so that 
61 = 8 2 = 0, we have no entanglement at all. In this case the ground state has energy -2 J and is doubly degenerate. 
For small, but equal, angles 6\ =6*2, there are two energy levels, one of energy value -2 J and the other with energy 
close to -2 J, with corresponding states, the first one equal to the Bell state (|01) — 1 10) ) j\[2 and the other close to 
the Bell state ( 1 01) + 1 10)) /V2. Thus, we get a maximally entangled qubit pair. Hence, we conclude that even a very 
small component of the magnetic field along the x-axis is adequate to create entangled states. In all cases we observe 
that the concurrence drops to zero for very strong fields. This is expected as for very strong fields the spins will be 
completely aligned along the field direction and hence the entanglement will drop to zero. It is worth noting that 
the smaller the x-component of the magnetic field becomes, the faster the concurrence drops to zero for B > 2 J. In 
the case where the magnetic field is along the x-axis the concurrence can be calculated analytically from the density 



matrix of the pure (ground) state, and we get C\ 2 = 1 + (B/J) 2 , see short dashed line {9\ = 9 2 = 7r/2) in Fig. 

1(a). The physical explanation of this behavior comes from the fact that the field tends to align the qubit spins in 
a different disentangled state from the spin-spin coupling. This implies that it is the trade off between the field and 
the Ising interaction that produces the entanglement. 

We now investigate the case of non-uniform magnetic field, which is actually more realistic than the uniform 
magnetic field case. In Fig. 1(b) we observe that for B 2 slightly different from B\ (here, B 2 = 1.0005[i?i = B}) the 
QPT disappears in cases of small x-component of the magnetic fields. Now the concurrence starts from vanishing value 
reach a maximum value for B < 2 J and then drops to zero very abruptly for B > 2 J. For small, but equal, angles 
9 1 =62, there are two energy levels one of energy value -2 J and the other with energy value close to -2Jbut in this 
case the corresponding states are not close to the Bell states ( 1 01) ± 1 10)) j\f2. As the x-component of the magnetic 
field increases the change from the vanishing value of the concurrence at zero field is very abrupt, as can be seen from 
the short dashed curve (9\ =62 = tt/2) in Fig. 1(b). For example, for 9\ = O2 = 0.3-7T and for a magnetic field of value 
B =0.01J the concurrence is practically unity. A very similar behavior is observed for slightly different directions of 
the magnetic fields i?iand B 2l as can be seen in Fig. 1(c). In this case, we observe that the QPT disappears for magnetic 
fields of the same magnitude but slightly different direction. It is important to point out that this phenomenon is 
present even for directions of the magnetic fields very close to the x-direction. Again, we get that there are two lowest 
energy levels corresponding to states that are not close to the Bell states ( 1 01) ± 1 10)) /y/2. For larger differences 
between the magnitudes of the magnetic fields for the two qubits, i.e. B 2 = 1.05i?i, the disappearance of the QPT 
is even more pronounced for cases characterized by small x-component of the magnetic fields [see solid curve in Fig. 
1(d)]. A systematic study of the field effects is depicted in Fig. 2, where we plot the concurrence for cases of small 
[Figs. 2(a) and (b)] and large [Figs. 2(c) and (d)] difference in the magnetic field magnitude. We observe that the 
entanglement achieved for low B x components is the one that it is more influenced. For example, even for a very small 
difference in the field amplitudes (of the order of 1.05) the maximum entanglement is diminished. For a larger value 
of the amplitudes of the field ratio, e.g. for 1.5, the concurrence practically disappears (not shown). The reason for 
this behavior is due to the fact that the eigenvectors are in this case the non-entangled pure states 1 01) , 1 10) and 1 00) . 
As we see from the same figure, the influence of the non-uniformity of the magnetic fields becomes smaller the more 
the direction of the fields get closer to a direction along the x-axis. In the same figure we investigate the influence of 
the difference in the field directions. We note that the maximum change in the behavior is obtained for the case that 
B x is small. 

We will now study the case of finite temperature, i.e. the case of thermal entanglement. We examine how the 
effects discussed in Fig. 1(a) changes for two different values of temperature with B\ = B 2 (= B), 9\ = 9 2 . From 
Figs. 3 and 4 we clearly see that the critical temperature depends on the orientation of the fields and that the field 
at which the maximum entanglement occurs as the temperature increases shifts to larger values of the magnetic field 
amplitudes. First, we study the case of the low B x component, where we observe a sharp peak in the concurrence. We 
note there is a fast drop of the concurrence to zero for a low B x component and for low values of the field magnitude. 
The behavior found in Figs 3 and 4 has already reported and discussed in Ref. [14]. From the same plots we observe 
that as the direction of the fields gets close to an orthogonal direction with respect to the Ising axis, the maximum 
entanglement occurs at lower values of the magnetic filed magnitude. In the case of 9 = tt/2, the concurrence does 
not drop to zero very fast as the eigenvectors of all eigenstates are either Bell states or a linear combination of Bell 
states. 
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Comparing Fig. 1(a) and Fig. 3(a), we clearly see that the QPT is not present in the finite temperature case. 
Moreover, we perform calculations for finite temperature with magnetic fields the same as those in Figs 1 (b) and 1 (c) . 
We have found that for these values of the external parameters, the concurrence is practically the same as in Fig. 3(a). 
This finding is valid once the magnitudes of the two applied fields are very similar. For case of larger difference in the 
field magnitude, the finite temperature results are more influenced as one sees by comparing Fig. 1(d) and Fig. 3(b). 
In this case of finite temperature under higher anisotropy of fields the QPT disappears, independently of the direction 
of the applied fields. Another, general feature seen from Fig. 4 (in comparison to Fig. 2) is that in all the cases the 
finite temperature leads to a decrease in the concurrence and shifts the maximum value of the concurrence at stronger 
applied fields. The effect is more pronounced for directions of the fields along an axis perpendicular to the Ising axis. 

It is known [14] that there is an angle Q\ = 62 = 9 (B, T) where the entanglement is maximum for a given temperature 
and amplitude (B\ = 82). This feature, known as the phenomenon of magnetically induced entanglement, has been 
explained heuristically assuming that with B x and B z fixed, the entanglement should change continuously with 
temperature. As the increase in the temperature widens the low-entanglement zone around B x = [14] and the 
entanglement has to fall for large B x , it is expected that at some intermediate value of B x , the maximal entanglement 
will be reached. The preferred angle traverses from — at zero temperature to 9 = ir/2 at T ~ J/ks- First, we 
investigate the zero temperature case. For parallel directions of the magnetic fields, the maximum concurrence is 
achieved for fields pointing along the x-axis (not shown here). Moreover, we observe that once we fix the direction 
of one of the applied fields the direction of the field applied at the other spin depends on the relative magnitude of 
the fields. This dependence is minimal once one of the fields is along the x-axis (not shown here). We note that 
there is a cutoff value of the magnetic field magnitude applied on the second spin, above which the relative angle 
remains practically the same. This cutoff value depends on the direction of the magnetic field applied at the first of 
the spins. We have shown that for non-uniform fields (B\ = J, B2 = 1.5 J) the maximum concurrence is obtained for 
non-parallel directions of the two applied fields. We observe that the smaller the B x component of one of the fields 
becomes the more the system deviates from the parallel direction. Actually, we have found that for fields of equal 
magnitudes and smaller than 2 J, no matter what is the direction of the field applied at the first spin, the maximum 
concurrence is achieved for magnetic field at the second spin pointing at the same direction as in the first one. This 
is not true for fields of equal magnitudes but with values larger than 2 J. The deviation from parallel fields becomes 
more pronounced the more the direction of one of the fields gets closer to the Ising direction. It is well understood 
that if both fields get very large values the concurrence is zero. In our study we found that there is no entanglement 
even if either one of the fields is very large. How large should the field be in order to destroy entanglement depends 
on the orientations of the fields. A rule of thumb is that we need larger fields as the directions of the fields gets closer 
to a direction perpendicular to the Ising direction. 

In Fig. 5(a) we confirm that for fields of equal magnitude (B\ = B2), the maximum concurrence occurs for 6\ = 62- 
On the contrary, Fig. 5(b) shows that for non- uniform fields the maximum entanglement occurs for unequal angles 
(#i 7^ #2)- Actually, as the difference in the magnitude of the two fields increases, the maximum entanglement occurs 
for direction of the larger field along the x-axis independent of the direction of the other field. A very similar behavior 
was found for a finite temperature as shown in Fig. 6. 

IV. CONCLUSIONS 

In summary, in the present work, a systematic investigation was performed on the entanglement in a two-spin Ising 
model in a site-dependent magnetic field. One of the most interesting results is the finding that the QPT observed at 
zero temperature when a uniform magnetic field is applied, disappears with the introduction of a very small difference 
in the applied fields. This difference could be either a difference in magnitude or a difference in the direction of the 
magnetic fields. Moreover, we have found that for parallel fields with direction close to the z-axis (Ising direction), 
small differences in the magnetic field magnitudes result in very weak entanglement. On the contrary, a very large 
asymmetry in the amplitudes of the applied fields has small effects on the well-entangled states obtained for fields with 
large x-component. We have also studied the phenomenon of the magnetically induced entanglement and observe that 
for equal magnitudes of the external magnetic fields the maximum concurrence occurs for parallel directions of the 
fields. Once a non-uniformity is introduced the maximally entangled states obtained for non-parallel orientations of 
the fields. In addition, we have shown that the concurrence drops to zero even if only one of the fields gets very large 
values. Finally, we have derived an analytic result valid for ground state entanglement and thermal entanglement. The 
analytic result, which has been confirmed numerically, predicts that we get vanishing concurrence once the direction 
of one of the fields is along the Ising direction. 



APPENDIX 



In the Appendix we discuss the special case of one magnetic field parallel to z-azis and the other in any direction. 
We find analytic expressions that predict that the concurrence in this case is always zero. The Hamiltonian studied 
in this case is given by 



H = J(o\al + erf of) + B x u\ + B 2 cos6<t* + B 2 sin 6 (erf + erf) /2 



(3) 



where a a are the Pauli operators and J is the strength of the Ising interaction. Also, B\ and B 2 are the magnitudes 
of the external magnetic fields. We assume that the first spin feels a magnetic field along the ^-direction and the 
second field has an arbitrary direction, defined by the angle 9 between the field and the Ising direction. 
The eigenenergies and the eigenstates can be calculated analytically. The eigenenergies are: 



= -B x ± [4 J 2 + B%- 4JB 2 cos 6] 1/2 

E ± = B 1 ±[AJ 2 + B 2 2 +AJB 2 cos6} 1/2 , 

The corresponding normalized eigenvectors are 

\X±) = o± |00) + b± |01) , 
\Y±) = c±|10)+d±|ll) , 



where, 



a± = 



B 2 smO 



>/2, 



[4 J 2 + B\ - AJB 2 cos 9} 1/4 (4 J 2 + B\ - 4JB 2 cos6») i/z ± (B 2 cos0- 2 J) 



1/2 



1/2 



b± = ± 



c± = 



(4 J 2 + B\ - 4 JB 2 cos 0) 1/2 ± (B 2 cos 9 - 2 J) 



1/2 



V2 



[4J 2 + B 2 - 4JB 2 cost 



,1/4 



B 2 sin 6 



V2 



[4 J 2 + B 2 + 4 JB 2 cos 0] 1/4 (4 J 2 + B 2 + AJB 2 cos (9) v ' ± (B 2 cos + 2 J) 



1/2 



1/2 



d± = ±- 



(4 J 2 + B\ + 4 JB 2 cos (?) 1/2 ± (S 2 cos 6 + 2 J) 



1/2 



V2 [4J 2 + B 2 + 4JB 2 costf] 1/4 

Now the density matrix is estimated by the following expression 

p = Z-\e-< 3E - \X_) (X_\+e^ E + \X+) (X+\ + e^ E - \Y_) (Y. 
where the partition function is defined as 

Z = e-? E - + e-? E + + e~P E - + e~^ E + . 



~ 0E + \y+) (y+w, 



(4) 
(5) 



(6) 
(7) 



(8) 

(9) 
(10) 

(11) 

(12) 
(13) 



and = 1/ksT. Then, the spin-flipped density matrix operator R 12 in the regular basis representation has the 
following form 



R 



12 



ACQ 

D B 

B -C 

-D A 



(14) 



where, A = pupa - p 12 p 34 , B = p 22 p 3 3 - P12P34, C = P33P12 - P11P34 and D = p i4 p 12 - p 22 p 3 4- The eigenvalues, A's 
of the spin flip density matrix operator are easily estimated as 



Ax = A 2 = [(A + B) + ((A - B) 2 + ACD) 1/2 ]/2 , 
A 3 = A 4 = \{A + B)- {{A-Bf +ACD) 1,2 ]/2. 



(15) 
(16) 



Hence, the concurrence defined as C12 = max{v^Vi~ — yA2 — V^3 — \/AtjO}, is always zero no matter what are the 
values of A, B, C and D. Therefore, we conclude that independently of the magnitudes of the magnetic fields and 
independently of the direction of one of the fields, if one of the fields points along the Ising direction there is no 
entanglement in the system. 

It is worth mentioning that the above analysis is greatly simplified for the zero temperature case as only the ground 
state is populated and hence p = |^o) (^o|> where is the index for the ground state. In this case the only non- 
vanishing terms in the density matrix operator are the pn, P22 and p\2- Then, it is rather straightforward to show 
that the spin-flip density operator matrix is the zero 4x4 matrix. 
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FIG. 1: Plots of concurrence as a function of the magnetic field magnitude B(=Bi) in units of J at zero temperature, (a) 
Equal magnetic field magnitudes (Bi = B2 = B). Solid curve is for 61 — 82 = 0.0I7T, long and short dashed curves correspond 
to 61 = 62 — O.Itt and 8i = 62 = 0.5-7T, respectively, (b) The same as (a) but with unequal field amplitudes (i?2=1.0005i3i). 
(c) Equal magnetic field magnitudes (Bi = B2 = B). Solid curve is for 61 = O.OI71-, 62 = 0.0117T, long dashed curve is for 
8\ = 0.l7r,#2 = 0.117T and short dashed curve is for 61 = 0.57T, #2 = 0.517T. (d) The same as (a) but with unequal field 
amplitudes (-B2=1.05i?i) 




FIG. 2: Plots of concurrence as a function of the magnetic field magnitude B(=Bi) in units of J at zero temperature, (a) 
Unequal field amplitudes (i?2=1.0005-Bi). Solid curve is for Oi = 0.017T, long and short dashed curves correspond to 0\ = 0.17T 
and 0i = 0.57T, respectively. Also, 62 — Oi = 0.017T. (b) Same as (a) but with 62—61 = 0.17T. (c) Same as (a) but with 
B 2 =1.05Bi. (d) The same as (c) but with 62 - 61 = 0.17T. 




FIG. 3: Plots of concurrence as a function of the magnetic field magnitude B(=Bi) in units of J at non-zero temperature (T= 
O.OU/fcs)- (a) Equal magnetic field magnitudes (Bi — B2 = B). Solid curve is for 61 — 62 = 0.017T, long and short dashed 
curves correspond to 61 = 62 = O.lir and 61 = 62 = 0.57T, respectively, (b) Same as (a) but with unequal field amplitudes 
(B 2 = 1.05Bi). 




FIG. 4: The same as in Fig. 2 but with non-zero temperature T= O.lJ/fcs. 




FIG. 5: Contour plots of concurrence as a function of the directions of magnetic fields 6\ and 82 (polar angles between field 
and Ising z-axis) for zero temperature. In both plots B\=2.1J . In (a) B2 = B\ and in (b) i?2=3Bi. 
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FIG. 6: Contour plots of concurrence as a function of the directions of magnetic fields 9\ and #2- In both plots T= U/fcs and 
Bi=2.1J. In (a) B 2 = Bi and in (b) B 2 =3Bi. 



